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REFOCUSING  A  RADIO  TELESCOPE  TO  IMAGE  SOURCES 
WITHIN  THE  NEARFIELD  OF  THE  ANTENNA  ARRAY 


1  THREE  FRAUNHOFER  APPROXIMATIONS  FOR  THE  PROPAGATION  OF  FIELD 
AMPLITUDE 

1.1  Introduction 

The  Fraunhofer  approximation  is  one  of  the  most  useful  mathematical  tools  used  in  optics.  It  is 
widely  used  by  optical  engineers  to  model  systems  and  by  researchers  to  model  physical  phenomena. 
The  Fraunhofer  approximation  is  also  widely  used  by  physicists  in  connection  with  radiation  prob¬ 
lems,  and  by  electrical  engineers  to  deal  with  antennas  and  electromagnetic  wave  propagation.  How¬ 
ever,  it  is  not  so  well  known  that  the  Fraunhofer  approximation  as  usually  used  in  these  different  dis¬ 
ciplines  is  different  in  some  important  ways.  These  differences  and  the  resulting  capabilities  to  model 
phenomena  have  given  rise  to  their  different  use. 

In  Section  1  we  examine  applications  of  three  Fraunhofer  approximations  to  simple  field  theory, 
where  we  deal  with  the  propagation  of  a  (generally  complex)  scalar  field  amplitude  from  sources  to 
detectors.  In  Section  2  we  extend  the  study  to  applications  of  the  Fraunhofer  approximations  to 
second-order  coherence  theory. 

The  present  section  derives  these  three  forms  of  the  Fraunhofer  approximation  by  using  similar 
notation  to  simplify  comparison.  Then  each  of  the  three  kinds  of  Fraunhofer  approximations  are 
examined  and  compared  to  the  others  to  point  out  the  differences  in  the  conditions  for  validity, 
mathematical  structure,  and  unique  capabilities  in  application. 

1.2  Detector  in  the  Farfield  of  the  Source  Distribution 

First  we  will  briefly  review  the  Fraunhofer  approximation  as  it  usually  appears  in  radiation 
problems  in  mathematical  physics  [1,  Section  9.1].  Here  we  establish  a  specific  set  of  conditions  and 
notation  to  contrast  this  formulation  with  different,  but  similar  Fraunhofer  approximations  frequently 
used  in  optics  and  in  electrical  engineering. 

Consider  the  configuration  shown  in  Fig.  1  with  a  three-dimensional  source  distribution  local¬ 
ized  to  a  region  about  the  origin  and  radiating  an  electromagnetic  field  toward  a  detector  located  a 
distance  R  from  the  origin  in  any  direction.  We  assume  that  the  source  is  a  primary,  three- 
dimensional  source  distribution  so  that  a  typical  scalar,  monochromatic,  Cartesian  component  of  the 
electromagnetic  field  i/<x)  (with  suppressed  exp  (-iwt)  time  dependence)  is  related  to  the  source  dis¬ 
tribution  p(x)  by  the  scalar  Helmholtz  equation 

(V2  +  k2)m  =  -4rp(x)  (1) 
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three-dimensional  source  distribution  about  the  origin 


[2,  Eq.  (2.8)],  where  k  =  2 t/X  =  w/c.  The  solution  to  Eq.  (1)  that  represents  an  outgoing  wave  at 
infinity  is 

0 o 

^(*,)  =  555  p<*d  exP  I  ri  l  v  I  ri  I  d\  -  (2) 

—  oo 

where  xj  and  X!  represent  radius  vectors  from  the  origin  to  the  source  and  field  points  respectively 
and  f[  represents  the  vector  X[  —  x|,  as  shown  in  Fig.  1.  Following  a  well-known  procedure  we 
expand  |  r(  |  into  the  binomial  series 

I r,  |  =  V(x,  -  x j)2  +  (y,  +  yj)2  +  (z i  -  zj)2 


~R  -  (x,  ■  x])//?  +  0(  |  xj  1 2/R), 


(3) 


where  R  =  |  X)  | .  If  the  source  distribution  is  confined  to  a  spherical  region  within  a  radius  a  of  the 
origin  and  the  detectors  are  all  far  outside  the  Rayleigh  range  of  the  source  so  that 

R  >>  ka2,  (4) 

then  we  may  approximate  the  value  of  |  rj  |  in  the  denominator  of  the  kernel  of  the  integral  in  Eq.  (2) 
by  the  first  term  in  Eq.  (3)  and  the  value  of  |  rj  |  in  the  phase  by  the  first  two  terms  to  obtain  the 
Fraunhofer  approximation  of  the  first  kind  given  by 


=  exp  (ikR)/R  JJJ  p(x|)  exp  (-iks  ■  x|)d3xj,  (5) 
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where  s  is  the  unit  vector  from  the  origin  to  the  field  point  at  x,  and  x{  is  the  three-dimensional  radius 
vector  from  the  origin  to  a  point  in  the  source  distribution. 

Equation  (5)  gives  the  Fraunhofer  approximation  for  the  radiation  field  from  a  primary,  three- 
dimensional  source  distribution  about  the  origin  as  it  often  appears  in  physics  problems  [1,  Eq.  (9.8)]. 
In  optical  diffraction  problems  a  similar  but  somewhat  different  form  of  Fraunhofer  approximation  is 
usually  used  [3,  Section  8.3.3].  Instead  of  a  primary,  three-dimensional  source  distribution  a  planar, 
secondary  source  distribution  is  usually  assumed  and  the  field  is  limited  to  the  z  2  0  half  space.  In 
such  problems  we  assume  that  the  primary  sources  are  distributed  over  the  z  <  0  half  space  and  radi¬ 
ate  toward  the  z  ~  0  plane.  We  assume  that  owing  either  to  the  nature  of  the  source  distribution  or 
to  an  aperture  in  the  z  =  0  plane  the  field  in  this  plane  vanishes  outside  of  some  domain  within  a 
radius  b  of  the  origin,  as  shown  in  Fig.  2.  Then  this  boundary  condition  for  the  field  over  the  z  =  0 
plane  ^(0)(xj)  is  treated  as  a  secondary  source  distribution. 


Fig.  2  —  Configuration  for  the  Fraunhofer  approximation 
of  the  second  kind  with  the  detector  in  the  farfield  in  the 
z  a  0  half  space  of  a  two-dimensional  source  distribu¬ 
tion  over  the  z  -  0  plane  near  the  origin 


For  a  secondary  source  distribution  the  homogeneous  Helmholtz  equation  describes  the  field 
over  the  z  >  0  half  space  [4,  Section  3.2].  A  well-known  solution  [4,  Section  3.4]  of  this  equation 
for  the  field  over  the  detectors  is  given  by  the  Rayleigh  integral 

ao 

Mx2)  =  - 1  /2t  jj  t<0>W)d/dz2[exp  (ik  |  r2  ] )/ 1  r2 1  )d2x2,  (6) 


where  x2  and  x2  are  the  radius  vectors  from  the  origin  to  the  field  point  at  the  detectors  and  the 
source  point  in  the  z  =0  plane,  respectively,  and  r2  is  the  vector  x2  -  x2,  as  shown  in  Fig.  2.  If  we 
assume  that  k  |  r2 1  is  sufficiently  large,  then  we  can  make  the  approximation 

-3/dz2[exp  (ik  |  r2 1 )/ 1  r2 1  ]  =  ( ~z2/ 1  r2 1  )(ik  -  1  / 1  r2 1  )[exp  (ik  |  r2 1 )/ 1  r2 1  ] 

~  ( ~ikz2/ 1  r2 1  )[exp  (ik  |  r2 1 )/ 1  r2 1  ].  (7) 
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It  is  useful  to  expand 

I r2 1  =  V(x2  -  x'2f  +  (y2  ~  yif  +  z$ 


~R  -  (x2  •  x{)/R  +0(1x11 2/R),  (8) 

where  R  has  the  same  meaning  as  in  Eq.  (3).  If  we  also  assume  that  the  secondary  source  in  the 
z  =  0  plane  is  confined  to  the  interior  of  a  circle  of  radius  b  around  the  origin  and  that  the  detectors 
are  all  located  outside  of  the  Rayleigh  range  of  the  source  distribution  so  that 

R  >>  kb2,  (9) 

then  upon  substitution  from  Eq.  (7)  into  Eq.  (6)  and  replacing  the  value  of  |  r2 1  in  the  denominator 
of  the  kernel  with  the  first  term  in  Eq.  (8)  and  the  value  in  the  phase  by  the  first  two  terms  in  Eq.  (8) 
we  have  the  Fraunhofer  approximation  of  the  second  kind  given  by 


\J/(Rs)  =  (-ik/2ir)  cos  0  (exp  ( ikR)/R )  fj  i/'(0)(x2)(exp  (-iks  ■  x2>f2x2,  (10) 


where  x2  is  a  vector  within  the  z  =  0  plane,  i.e.,  x2  =  (x2,  y2,  0),  s  is  the  three-dimensional  unit 
vector  from  the  origin  in  the  direction  of  x2,  and  0  is  the  angle  that  s  makes  with  the  +z  axis. 

The  Fraunhofer  formula  Eq.  (10),  for  the  field  owing  to  a  secondary  source,  is  frequently 
employed  in  optics  to  deal  with  diffraction  problems  [4,  Eq.  (4.13)].  The  relation  given  in  Eq.  (5) 
for  a  three-dimensional,  primary  source  differs  from  Eq.  (10)  more  than  one  might  suppose  in  view 
of  the  similarity  of  their  derivations.  The  main  mathematical  differences,  which  will  be  examined  in 
more  detail  in  the  next  subsection,  arise  from  the  fact  that  the  integral  appearing  in  Eq.  (5)  is  three- 
dimensional  unlike  the  integral  that  appears  in  Eq.  (10)  and  that  s  is  a  unit  vector  so  that  the  three- 
dimensional  Fourier  transform  appearing  in  Eq.  (5)  is  evaluated  only  on  a  sphere  of  unit  radius  in 
(sx ,  sy,  sz)  space.  It  is  well  known  that  these  mathematical  properties  give  rise  to  some  very  impor¬ 
tant  properties  of  radiating  source  distributions,  and  are  of  particular  interest  when  dealing  with 
inverse  problems.  This  subject  will  be  discussed  further  in  Subsection  1.4. 

A  third  form  of  the  Fraunhofer  approximation  exists  that  is  frequently  used  by  electrical 
engineers  in  dealing  with  antenna  problems.  It  is  also  closely  related  to  the  Fraunhofer  approxima¬ 
tions  appearing  in  Eqs.  (5)  and  (10)  but  employs  a  significantly  different  set  of  assumptions  and, 
therefore,  has  very  different  applications.  We  present  a  brief  derivation  of  this  Fraunhofer  approxi¬ 
mation  to  establish  a  precise  set  of  conditions  and  notation  before  we  consider  these  different  formula¬ 
tions  in  more  detail. 

1.3  Source  Distribution  in  the  Farfield  of  the  Detector  Array 

Figure  3  shows  a  configuration  with  a  distribution  of  detectors  located  near  the  origin  in  the 
z  =  0  plane  that  receive  electromagnetic  radiation  from  sources  located  a  distance  R  from  the  origin 
in  the  z  2:  0  half  space.  Consider  only  a  secondary  source  distribution  (i.e.,  a  boundary  condition 
for  the  field  over  the  hemisphere)  rather  than  a  true  primary  source  distribution.  Comparison  of  Figs. 
2  and  3  shows  that  the  two  configurations  are  identical  except  for  interchanging  the  positions  of  the 
sources  and  the  detectors.  Because  of  this  similarity  the  fields  at  the  sources  and  detectors  can  be 
related  by  using  Eq.  (10)  and  a  well-known  reciprocity  theorem  from  antenna  theory. 
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Fig.  3  —  Configuration  for  the  Fraunhofer  approximation 
of  the  third  kind  with  the  sources  outside  of  the  Rayleigh 
range  in  the  z  2  0  half  space  of  a  two-dimensional 
detector  array  over  the  z  =  0  plane  near  the  origin 


This  reciprocity  theorem  [5,  Section  10.2]  leads  to  the  well-known  conclusion  that  the  transmit¬ 
ting  and  receiving  patterns  of  any  antenna  are  the  same.  This  theorem  can  be  applied  directly  to 
transmitting  patterns  such  as  that  given  by  Eq.  (10)  to  give  the  equivalent  receiving  pattern  of  a  detec¬ 
tor  array  over  the  z  =  0  plane.  Thus,  if  we  assume  that  the  detectors  are  all  confined  to  a  region  of 
the  z  =  0  plane  within  a  circle  of  radius  b'  from  the  origin  and  that  the  sources  of  the  field  all  are 
outside  of  the  Rayleigh  range  of  the  source  distribution,  i.e., 

R  >>  kb'2,  (11) 

then  Eq.  (10)  can  be  reinterpreted  as  an  inverse  propagation  integral  (after  taking  the  complex  conju¬ 
gate  of  the  kernel  as  discussed  in  Appendix  A).  Thus  we  obtain  the  field  at  a  particular  source  point 
i£(x3)  in  the  farfield  of  the  detector  as  a  function  of  field  distribution  i(0\x^)  over  the  detectors  in  the 
z  =  0  plane. 

It  is  well  known  that  the  Fraunhofer  approximation  of  the  second  kind  cannot  simply  be  inverted 
to  completely  determine  the  source  distribution  over  the  z  =  0  plane  from  measurements  of  the  far¬ 
field.  As  discussed  in  detail  in  Subsection  1.4,  this  results  from  the  fact  that  s  is  a  unit  vector  that  is 
only  defined  in  Eq.  (10)  for  directions  from  the  origin  into  the  z  S  0  half  space.  However,  this 
inversion  problem  disappears  when  Eq.  (10)  is  reinterpreted  by  use  of  the  reciprocity  theorem  to  be  a 
receiving,  rather  than  a  transmitting,  antenna  pattern  and  the  sources  of  the  field  are  defined  to  be 
located  only  in  directions  from  the  origin  to  the  z  >  0  half  space,  instead  of  on  the  z  =0  plane. 
Since  the  integral  in  Eq.  (10)  is  a  Fourier  transform  and  since  the  function  ^{R s)  is  defined  to  be 
nonvanishing  only  for  s  into  the  z  >  0  half  space  (see  Subsection  1.4),  Eq.  (10)  can  be  easily 
inverted  to  give  the  Fraunhofer  approximation  of  the  third  kind 

^<0)(Xj)  =  i-ik/lr)  Jf  [4iR$)exp(ikR)/R]exp(-iks  - x3)R2(dsxdsy/sz)  (12) 
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where  s  =  (sx,  sy,  sz)  is  the  radius  vector  from  the  origin  to  the  source  point  at  x3  and  x3  is  the 
radius  vector  from  the  origin  to  a  point  within  the  detector  array.  In  Eq.  (12)  R  is  generally  a  func¬ 
tion  of  s  since  we  have  had  to  make  no  mathematical  restriction  that  the  surface  containing  the 

sources  be  exacdy  hemispherical.  Thus  R(s)  cannot  be  taken  outside  of  the  integrals  in  Eq.  (12). 

•  «  > 

The  Fraunhofer  formulas  (Eqs.  (5),  (10),  and  (12))  do  not  have  the  same  conditions  for  validity. 
The  Fraunhofer  approximations  given  by  Eqs.  (5)  and  (10)  require  that  the  source  distribution  be 
localized  about  the  origin  of  coordinates  and  that  the  detectors  all  be  in  the  farfield  of  the  source  dis¬ 
tribution.  For  Eq.  (5)  the  detector  may  be  in  any  direction  from  the  source;  however,  in  Eq.  (10)  the 
detector  must  be  located  in  the  z  ^  0  half  space.  The  different  Fraunhofer  formulation  given  by  Eq. 
(12)  requires  that  the  detectors  all  be  localized  about  the  origin  of  coordinates  and  that  the  sources  all 
be  in  the  farfield  of  the  detector  array  within  the  z  ^  0  half  space.  Thus  the  conditions  for  Eq.  (10) 
are  reciprocal  to  those  for  Eq.  (12).  Sometimes  both  the  source  distribution  and  the  detector  array 
are  localized  and  are  each  in  the  farfield  of  the  other  so  that  these  formulas  differ  only  by  a  coordi¬ 
nate  transformation  from  the  origin  within  the  source  distribution  to  the  origin  within  the  detector 
array.  However,  this  is  generally  not  the  case.  Usually  only  one  set  of  conditions  is  met  so  that  only 
one  of  the  Fraunhofer  approximations  is  valid. 

1.4  Differences  Between  Fraunhofer  Approximations 

As  we  noted,  the  three  kinds  of  Fraunhofer  approximations  that  we  derived  are  different  in  some 
significant  respects.  We  now  consider  them  one  at  a  time  and  point  out  some  of  their  unique 
features. 

Only  the  Fraunhofer  approximation  of  the  first  kind  in  Eq.  (5),  involves  a  three-dimensional 
Fourier  transform,  which  has  some  interesting  consequences.  Although  the  transform  is  three  dimen¬ 
sional,  the  domain  of  ^(Rs)  in  Eq.  (5)  is  only  two  dimensional.  Since  s  =  (sx,  sy,  sz)  is  a  unit  vec¬ 
tor,  its  components  are  related  by 


sx  +  Jy2  +  sz2  =  1.  (13) 

Thus  the  Fourier  transform  in  Eq.  (5)  has  clear  physical  meaning  only  for  values  of  (sx,  sy,  sz  )  over 
the  surface  of  a  sphere  in  s  space.  As  a  consequence,  the  farfield  radiation  from  a  three-dimensional 
source  distribution  is  completely  specified  by  only  the  spatial  Fourier  components  of  the  source  that 
have  a  period  of  one  wavelength  [6,  Theorem  !]•  This  property  has  proven  useful  in  studying  the 
question  of  nonradiating  source  distributions  16-10],  scattering  problems  [11],  and  inverse  scattering 
[12].  It  appears  likely  that  similar  formulations  could  also  be  applied  to  phase  matching  problems. 
These  physical  properties  of  three-dimensional  source  distributions  are  not  described  by  the 
Fraunhofer  approximations  of  the  second  or  third  kinds. 

Another  interesting  consequence  of  the  three-dimensional  integral  in  Eq.  (5)  involves  an  impor¬ 
tant  class  of  inverse  problems.  These  inverse  problems  employ  measurements  of  farfield  data  to  calcu¬ 
late  the  source  distribution.  We  noted  from  Eq.  (5)  that  farfield  measurements  give  the  Fourier 
transform  of  the  three-dimensional  source  distribution  only  over  the  sphere  of  unit  radius  from  the 
origin  in  a  three-dimensional  s  space.  To  determine  the  source  distribution  for  an  arbitrary  source  by 
Fourier  inversion  of  this  data,  the  data  would  have  to  be  known  everywhere  over  s  space.  Some  spe¬ 
cial  cases  are  known  for  which  farfield  measurements  are  sufficient  to  reconstruct  the  source  distribu¬ 
tions.  In  scattering  problems  using  the  first  Bom  approximation,  for  example,  the  source  distribution 
represents  a  product  of  the  known  incident  field  and  the  scattering  potential  of  the  scatterers.  There¬ 
fore,  the  scattering  potential  can  be  computed  from  a  series  of  farfield  measurements  with  different 
incident  fields  [13,  14], 
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The  Fraunhofer  approximation  of  the  second  kind  given  by  Eq.  (10)  is  different  in  two  obvious 
ways  from  Eq.  (S).  Equation  (10)  contains  constant  and  angle  factors  that  result  from  the  use  of  a 
secondary  source  distribution  and  a  two-dimensional  Fourier  transform. 

In  Eq.  (10)  the  Fourier  transform  is  from  the  source  plane  (x^,  y 2)  to  the  plane  (sx,  sy),  where 
the  third  component  of  the  unit  vector  s  is  given  by 

sz  =  cos  8  =  V 1  —  sx  —  Sy,  if  s }  +jy  s  1.  (14) 

Fourier  components  of  ^(Xj)  inside  of  the  unit  circle  in  Eq.  (14)  for  which  sz  is  defined  are  called 

low-frequency  spatial  frequency  components  of  the  source  and  those  outside  are  called  high  spatial 
frequency  components. 

This  difference  in  dimensionality  gives  rise  to  the  description  of  some  interesting  physical 
phenomena  that  are  not  described  by  the  Fraunhofer  approximation  of  the  first  kind  in  Eq.  (5).  An 
arbitrary  source  distribution  contains  high  spatial  frequency  components  outside  of  the  unit  radius  in 
(sx,  sy)  space  for  which  sz  is  defined  in  Eq.  (14).  These  high-frequency  spatial  frequency  com¬ 
ponents  have  no  effect  on  the  farfield  in  any  direction  from  the  origin  into  the  z  >  0  half  space. 

They  are  associated  with  evanescent  plane  waves  that  do  not  propagate  into  the  farfield.  It  is  possible 
to  reconstruct  the  source  distribution  by  using  only  the  low-frequency  spatial  frequency  components 
that  are  available  from  farfield  measurements.  However,  the  use  of  only  low-frequency  components 
results  in  a  spatial  resolution  limit  similar  to  that  imposed  by  a  diffraction  limited  lens  ( NA  =  1)  to 
detail  within  the  source  distribution  separated  by  intervals  larger  than  a  wavelength. 

The  Fraunhofer  approximation  of  the  third  kind  in  Eq.  (12)  is  different  from  that  of  the  second 
kind  in  Eq.  (10)  in  that  the  Fourier  transform  is  from  the  s  space  to  xj  space  instead  of  the  other  way 
around.  In  Eq.  (12)  it  is  the  sources  that  lie  in  directions  away  from  the  origin  in  the  z  SO  half 
space  and  the  detectors  of  the  field  that  are  in  the  z  =0  plane.  The  range  of  integration  in  Eq.  (12) 
is  only  over  real  angles  into  the  z  s  0  half  space  so  that  the  x  and  y  components  of  s  are  bounded  to 
the  interior  of  the  unit  circle,  as  given  by  Eq.  (14).  Thus  from  Eq.  (12)  we  see  that  the  field  over  the 
z  —  0  plane  contains  no  spatial  frequency  components  with  periods  smaller  than  a  wavelength. 
Furthermore,  by  using  a  well-known  theorem  of  Fourier  analysis  we  know  that  the  field  over  the 
2=0  plane  can  be  extended  as  an  entire  function  of  complex  x$.  Since  the  secondary  source  distri¬ 
bution  is  only  defined  physically  for  x  and  y  components  of  s  within  the  unit  circle  we  can  define  the 
function  to  be  zero  elsewhere.  With  this  definition  there  is  no  problem  in  inverting  the  Fourier 
transform  in  Eq.  (10),  which  is  interpreted  as  a  receiving  antenna  pattern,  to  obtain  Eq.  (12).  A 
consequence  of  these  properties  of  the  Fourier  transform  in  Eq.  (12)  is  that  there  appears  to  be  no 
limitation  in  the  angular  resolution  to  which  the  source  distribution  can  be  calculated  other  than  that 
arising  from  the  limitation  on  the  size  of  the  detector  array  in  the  z  =0  plane. 

Both  the  Fraunhofer  approximations  of  the  first  and  second  kinds  are  valid  only  if  the  detectors 
are  in  the  farfield  of  the  source,  i.e.,  the  source  distribution  must  be  bounded  so  that  it  has  a  well- 
defined  Rayleigh  range,  and  the  detector  must  be  well  outside  of  this  Rayleigh  range.  In  many 
important  applications  the  source  distribution  is  either  not  bounded  or  too  large  for  this  approxima¬ 
tion.  A  radio  antenna  is  an  example.  Here,  the  domain  for  the  source  distribution  might  completely 
surround  the  antenna.  However,  if  the  detector,  like  a  radio  antenna,  is  bounded,  then  it  has  a  Ray¬ 
leigh  range  itself  and  for  sources  outside  of  this  Rayleigh  range  the  Fraunhofer  approximation  of  the 
third  kind  holds. 

An  important  difference  in  the  application  of  the  Fraunhofer  approximations  of  the  second  and 
third  kinds,  which  both  involve  calculating  the  field  from  a  known  planar  source  distribution,  can  be 
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shown  by  use  of  an  example.  Consider  a  planar  source  distribution  localized  within  a  circle  of  radius 
b  about  an  origin  at  O  in  the  plane  P  in  Fig.  4.  Assume  that  it  is  radiating  toward  a  planar  detector 
array  localized  withm  ?  radius  b '  of  an  origin  O '  in  the  plane  P '  a  distance  D  from  P . 


1 

*  V 

I 

b 

- i _ j 

* 

b' 

4 

0 

— ■  ■  _ _ _  r\ 

O' 

u - - - * 

Fig.  4  —  Notation  for  the  example  contrasting  calculations  per¬ 
formed  using  the  Fraunhofer  approximations  of  the  second  and 
third  kinds 


If  b  and  R  satisfy  Eq.  (9)  so  that  Eq.  (10)  is  a  valid  approximation,  then  to  calculate  the  field 
over  the  detector  array  from  a  known  source  distribution,  we  perform  a  two-dimensional  Fourier 
transform  of  the  source  distribution  data,  multiply  the  data  by  the  constants  and  s-dependent  factors 

outside  of  the  integral  in  Eq.  (10),  and  finally  transform  from  angular  to  linear  coordinates  1'  within 
the  P'  plane,  i.e.. 


i;  =  D(sx/yl  1  -  s?  -  jv2) 

=  D(Sy/yj  1  -  S?  -  Sy).  (15) 


Figure  5(a)  summarizes  these  steps. 

If  instead  b’  and  R  satisfy  Eq.  (11)  so  that  Eq.  (12)  is  a  valid  approximation,  then  to  calculate 
the  field  over  the  plane  P '  from  a  known  source  distribution  we  must  use  a  different  procedure.  First 

we  transform  the  source  distribution  from  linear  coordinates  |  in  the  P  plane  into  the  angular  coordi¬ 
nates  from  the  origin  at  O',  i.e.. 


*  =  {,/Vd2  +  g  +  $2 

and 

=  £y/VD2  +  $2  +  £2.  (16) 

Next  we  multiply  the  data  by  the  s-dependent  factors,  and  finally  we  Fourier  transform  the  resulting 
data  over  the  (sx ,  sy)  plane  to  obtain  the  field  over  the  £'  plane.  Figure  5(b)  shows  this  procedure. 
Study  of  this  figure  readily  shows  the  differences  in  the  two  procedures. 

Finally  we  consider  the  interesting  case  for  which  b,  b',  and  R  satisfy  both  Eqs.  (9)  and  (1 1)  so 
that  Eqs.  (10)  and  (12)  are  both  valid  approximations.  For  this  case  either  processing  scheme  shown 
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(a)  Field  computation  if  detectors  are  in  the  farfield  of  the  source 


(b)  Field  computation  if  sources  are  in  the  farfield  of  the  detector  array 


Fig.  5  —  Differences  in  propagation  calculations  using  the  Fraunhofer 
approximations  of  the  (a)  second  and  (b)  third  kinds 


in  Fig.  5  should  work.  Examination  of  the  nonlinear  coordinate  transformations  given  by  Eqs.  (IS) 
and  (16)  indicate  that  in  the  limit  as  R  gets  very  large  relative  to  both  b  and  b '  we  get 

tx  ~  t>sx, 


-  Dsy,  (17) 

and 

>X  ~ 

sy  -  i;/D,  (18) 

so  that  in  this  limit  the  transformations  represent  only  a  linear  scaling.  Substitution  from  Eqs.  (17) 
and  (18)  into  Eqs.  (10)  and  (12)  while  setting  R  =  D  gives 

oo 

\Kl)  =  -i*/2r[exp  (ikD)/D]  ff  ^(?)[exp  (-it?  •  f/D)d^xd^],  (19) 

—  OO 


for  either  case  which  shows  that  they  are  identical  in  this  limit. 

The  Fraunhofer  approximations  of  the  second  and  third  kinds  have  another  interesting  difference 
when  they  are  applied  to  inverse  problems.  As  we  discussed,  the  Fraunhofer  approximation  of  the 
second  kind  can  be  inverted  for  use  in  computing  the  source  distribution  from  measured  farfield  data, 
but  only  to  within  the  resolution  limit  of  a  diffraction-limited  lens  with  a  numerical  aperture  of  unity. 
The  Fraunhofer  approximation  of  the  third  kind  can  also  be  inverted  and  applied  to  inverse  problems, 
but  because  of  the  requirement  that  the  detector  array  be  limited  to  circle  of  radius  where  Eq.  (11)  is 
satisfied,  the  field  cannot  be  measured  outside  of  this  region.  This  restriction  leads  to  a  different 
resolution  limit  on  the  calculations  for  the  source  distribution.  In  this  case  the  fundamental  resolution 
limit  is  equivalent  to  that  of  a  diffraction-limited  lens  with  a  numerical  aperture  of  b'lR ,  which, 
because  of  Eq.  (11),  is  much  poorer  than  the  best  resolution  possible  when  using  the  Fraunhofer 
approximation  of  the  second  kind. 
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Table  1  summarizes  the  different  properties  of  the  three  kinds  of  Fraunhofer  approximations  dis¬ 
cussed  in  this  section. 


Table  1  —  Differences  Between  Fraunhofer  Approximations 


First  Kind 

Second  Kind 

Third  Kind 

Source  domain 

Near  origin 

z=0  plane 

2  >  0  half  space 

Source  dimen. 

3-D,  linear 

2-D,  linear 

2-D,  angular 

Det.  domain 

Far  from  origin 

Far  from  0 ,  z  >  0 

2=0  plane 

Det.  dimensions 

2-D,  angular 

2-D,  angular 

2-D,  linear 

FT  domain 

Sphere 

Plane 

Plane 

FT  dimensions 

3-D 

2-D 

2-D 

Large  parm. 

(Source  radius)  2/R 

(Aperture  dia.)  2/R 

(Det.  dia.)  2/R 

1.5  Conclusions 

The  Fraunhofer  approximation  can  take  at  least  three  forms  that  are  very  different  and  are  useful 
in  different  applications.  The  Fraunhofer  approximations  of  the  first,  second,  and  third  kinds,  defined 
in  this  section  are  well  known  in  radiation  physics,  optics,  and  electrical  engineering,  respectively. 
This  appears  to  be  because  their  different  properties  best  suit  them  for  important  applications  in  these 
fields.  However,  each  is  not  as  well  known  outside  of  its  respective  field.  Since  occasionally  prob¬ 
lems  arise  in  optics  that  can  be  best  analyzed  by  a  Fraunhofer  approximation  other  than  that  of  the 
second  kind,  it  would  seem  useful  to  be  aware  of  all  three  of  these  forms. 

2  THREE  FRAUNHOFER  APPROXIMATIONS  FOR  THE  PROPAGATION 
OF  THE  CROSS-SPECTRAL  DENSITY  FUNCTION 

2.1  Introduction 

Section  1  (also  Ref.  15)  shows  that  at  least  three  different  kinds  of  Fraunhofer  approximations  in 
scalar  diffraction  theory  have  important  differences  in  both  conditions  for  validity  and  mathematical 
structure. 

In  the  next  three  subsections  we  apply  each  of  these  approximations  to  coherence  theory  to 
describe  the  propagation  of  the  cross-spectral  density  function  from  the  sources  of  the  field  to  detec¬ 
tors  far  away.  As  was  the  case  in  diffraction  theory  we  find  significant  differences  between  the  three 
approximations.  In  the  fourth  subsection  we  show  that  in  two  of  these  approximations  the  propaga¬ 
tion  equations  can  be  easily  inverted  to  deal  with  inverse  problems.  The  inverse  propagation  equa¬ 
tions  suggest  a  method  for  refocusing  a  radio  telescope  to  observe  objects  in  the  nearfield  of  its 
antenna  array. 

2.2  Fraunhofer  Approximation  of  the  First  Kind 

Consider  Fig.  6  that  shows  a  three-dimensional  source  distribution  localized  to  a  region  about 
the  origin  and  radiating  an  electromagnetic  field  in  all  directions.  We  consider  the  typical  mono¬ 
chromatic,  Cartesian  components  of  the  field  ^(x,  w)  and  the  charge-current  distribution  p(x,  u), 
which  are  related  by  the  inhomogeneous  Helmholtz  equation 

(V2  +  k2)\f>(x,  w)  =  -4*p(x,  «),  (20) 

where  u  is  the  radial  frequency  and  k  =  2x/X  =  «/e . 
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Fig.  6  —  Configuration  for  the  Fraunhofer 
approximation  of  the  first  kind 


In  this  section  we  assume  that  the  source  fluctuates  in  an  irregular  manner  and  treat  u)  and 
p(x,  w)  (for  each  frequency  to)  as  random  variables.  The  second-order  statistical  properties  of  these 
random  variables  will  be  characterized  in  a  manner  conventional  in  coherence  theory  by  the  cross- 
spectral  density  functions  Wp  and  defined  by 

<p(xu  u)p*(x2,  to')  >  =  wp(xu  x2,  to)6(<o  -  to')  (21a) 

and 

<tf(x„  toW*(x2,  to')  >  =  W£\x,  x2,  to)5(to  -  to'),  (21b) 


which  represent  the  ensemble  average  of  the  random  variables  at  two  different  points  Xj  and  x2  and  at 
two  different  frequencies  to  and  to'. 

We  found  in  Subsection  1.2  that  if  the  source  is  completely  contained  within  a  sphere  of  radius 
a  centered  on  the  origin  and  if  the  field  point  is  at  a  distance  R  from  the  origin  that  is  outside  of  the 
Rayleigh  range  of  the  source,  i.e.. 


R  >>  ka2,  (22) 

then  from  Eq.  (S)  the  field  amplitude  at  x  is  given  as  a  function  of  the  three-dimensional  Fourier 
transform  of  the  source  distribution  over  the  source  points  x'  by  the  Fraunhofer  approximation  of  the 
first  kind,  viz. , 

OP 

'HRs)  =  [exp  (ikR)/R]  JJJ  p(x')  exp  (-iks  ■  x')d\\  (23) 

—  00 

where  R  =  |  x  |  is  the  distance  and  s  is  a  unit  vector  in  the  direction  from  the  origin  to  the  field  point 
at  x.  Upon  substitution  from  Eq.  (23)  into  Eq.  (21b)  and  then  by  using  Eq.  (21a)  we  find  that  the 
cross-spectral  density  function  in  the  farfield  of  the  source  is  given  by 

»V(*1  *1,  R2*2)  =  (CXP  [«*(*!  -  R2)]/R  1*2) 


x  fJJ  ill  *2)  e*P  !-«*(*!  •  *i  -  *2  •  xi)Jd3x| d*xi,  (24) 

—  00  —  00 
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where  we  drop  the  argument  u  from  here  on  (for  a  similar  result  see  Ref.  2,  Eq.  (3.3)).  Equation 
(24)  is  the  second-order  (in  the  field  amplitude)  form  of  the  Fraunhofer  approximation  of  the  first 
kind  that  gives  the  field  cross-spectral  density  function  in  the  radiation  zone  of  a  three-dimensional 
source  as  a  double  three-dimensional  Fourier  transform  of  the  cross-spectral  density  function  of  the 
source.  In  Eq.  (24)  s,  and  $2  represent  unit  vectors  from  the  origin  toward  the  field  points  at  x,  and 
x2,  respectively,  R{  =  |  x,  | ,  /?2  =  I  *2 1  >  the  two  volume  integrals  are  over  the  source  points  at 
xj  and  x2,  as  shown  in  Fig.  6. 

Let  us  now  consider  the  important  special  case  for  which  the  source  is  quasi-homogeneous  [16, 
Eq.  (3.5)]  so  that 


Wp(x[,  *$)  =  /p[(xj  +  x2)/2] /t(xj  -  x2),  (25) 

where  Ip  is  a  slow  function  relative  to  fi.  By  this  we  mean  that  n  is  zero  outside  of  some  domain  in 
which  |  xj  —  x2 1  is  smaller  than  some  fixed  constant  d '  and  Ip  is  constant  over  any  domain  where 
| xj  +  x2 1  /2  varies  by  no  more  than  d'  [16,  Section  3],  The  function  #t(x!_)  appearing  in  Eq.  (25)  is 
a  normalized  correlation  coefficient  that  has  a  maximum  value  of  unity  that  always  occurs  when 
|  x'_  |  =0.  By  substitution  from  Eq.  (25)  into  Eq.  (24)  we  obtain 

HV*lSl,  R&L>  =  feXP  (**(*!  -  ^2)]/^1^2ly"p(Sl  ~  ®2>£[(®1  +  S2>/2],  (26) 

where  the  caret  indicates  three-dimensional  Fourier  transforms 

ao 

/p(?)  =  JJJ  /p(x+)  exp  (-ikx'+  ■  £)d*x'+,  (27) 

-  OB 

and 

00 

A(I)  =  M*-)  exp  (-iJfcxl  •  ®d3xl  (28) 

—  OD 

[16,  Eq.  (2.3.10)]. 

Equation  (26)  is  an  important  result.  It  applies  to  the  usual  thermal  sources  and  to  most  other 
nonlaser  sources.  It  leads  to  a  pair  of  reciprocity  theorems  that  are  very  important  in  radiation  theory 
(see  Ref.  16,  Section  3  for  a  detailed  discussion  of  these  theorems). 

2.3  Fraunhofer  Approximation  of  the  Second  Kind 

Assume  that  the  source  is  a  planar,  secondary  source  in  the  z  =  0  plane,  as  shown  in  Fig.  7, 
instead  of  the  three-dimensional,  primary  source  distribution  discussed  in  the  last  subsection,  and 
assume  that  it  radiates  only  into  the  z  2:  0  half  space.  We  found  in  Subsection  1.2  that  if  such  a 
source  is  completely  contained  within  a  circle  of  radius  b  about  the  origin  and  if  the  field  point  is  at  a 
distance  R  from  the  origin  that  is  outside  of  the  Rayleigh  range  of  the  source,  i.e., 

R  >>  kb2,  (29) 

then  from  Eq.  (10)  the  field  in  any  direction  s  from  the  origin  into  the  z  £  0  half  space  can  be  given 
as  a  function  of  the  two-dimensional  Fourier  transform  over  the  source  distribution  by  the  Fraunhofer 
approximation  of  the  second  kind,  viz.. 


W?s)  =  (-ik/lx)  cos  6  exp  (ikR)/R  |J  ^(0)(x')  exp  (—/As  •  x'ydhi'. 


(30) 
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Fig.  7  —  Configuration  for  the  Fraunhofer 
approximation  of  the  second  kind 


where  $  is  the  angle  that  x  makes  with  the  +z  axis,  R  =  |  x  | ,  and  s  =  x/R .  By  substitution  from 
Eq.  (30)  into  Eq.  (21b)  and  then  by  using  Eq.  (21b)  we  have 

»V(J?,s„  ^2*2)  =  (k/2r)2  cos  0,  cos  02{exp  [t)k (Rx  -  R2)]/RlR2i 


00  00 

x  JJ  JJ  W'f  (x'„  xi)  exp  [~ik{ s,  •  xj  -  Sj  •  x^Jd^jd^.  (31) 

—  oo  —  oo 

Equation  (31)  gives  the  field  cross-spectral  density  function  in  the  radiation  zone  of  the  two- 
dimensional  source  distribution  as  a  double  two-dimensional  Fourier  transform  of  the  cross-spectral 
density  function  of  the  field  between  the  source  positions  xj  and  x2  in  the  *  =0  plane.  This  is  the 
second-order  form  of  the  Fraunhofer  approximation  of  the  second  kind. 

Now  we  will  consider  the  special  case  for  which  the  two-dimensional  source  distribution  is 
quasi-homogeneous,  i.e.. 


x})  =  /[(xj  +  x2)/2]  h(x\  -  xj),  (32) 

where  again  /  is  a  slow  function  relative  to  p  and  n  is  a  correlation  coefficient.  Upon  substitution 
from  Eq.  (32)  into  Eq.  (31)  we  have 

^(/?lS„  /JjSj)  =  (k  Hr)2  cos  0,  cos  02[exp  ik  (f?,  -  R2)/RXR2\ 


x  /( s,  -  s^jtKs,  +  S2)/2],  (33) 

where  the  tilde  represents  a  two-dimensional  Fourier  transform 

00 

/(I)  =  JJ  I(x'+)  exp  (~ikx'+  •  1)dzx'J,,  (34) 

—  00 

and 


00 

j*(?)  =  JJ  /»(*'-)  exp  (-ikxL  ■  fidhiL. 


(35) 
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Equation  (33)  is  very  important  in  coherence  theory.  It  applies  to  thermal  as  well  as  most  other 
nonlaser  sources.  It  leads  to  the  van  Cittert  Zemike  theorem  as  well  as  to  a  closely  related,  recipro¬ 
cal  theorem.  This  has  been  discussed  in  detail  elsewhere  [17,  Section  IV]. 

Equation  (31)  appears  to  be  mathematically  similar  in  form  to  Eq.  (24)  except  for  two  differ¬ 
ences.  The  Erst  difference  is  that  the  cross-spectral  density  function  inside  of  the  integral  in  Eq.  (31) 
is  between  values  of  the  field  amplitude  in  the  source  plane  rather  than  values  of  the  source  distribu¬ 
tion  itself.  This  gives  rise  to  some  additional  (obliquity)  factors  appearing  in  Eq.  (31).  The  second 
difference  is  that  the  Fourier  transform  in  Eq.  (31)  is  two  dimensional  rather  than  three  dimensional. 
As  we  shall  see  later,  this  difference  in  dimensionality  has  important  consequences. 

The  conditions  for  validity  of  Eqs.  (24)  and  (31)  are  of  course  quite  different.  Equation  (24) 
applies  to  the  Held  outside  of  the  Rayleigh  range  in  any  direction  from  a  three-dimensional  primary 
source  distribution  about  the  origin.  Equation  (31)  applies  to  the  field  outside  of  the  Rayleigh  range 
in  the  z  >0  half  space  of  a  secondary  planar  source  in  the  z  =  0  plane. 

2.4  Fraunhofer  Approximation  of  the  Third  Kind 

For  the  last  case,  we  assume  that  the  detector  array,  not  the  source,  is  distributed  over  the 
z  =  0  plane  near  the  origin  as  shown  in  Fig.  8.  Then  we  assume  that  it  is  the  source  of  the  field  that 
is  located  in  the  z  20  half  space  away  from  the  origin.  Comparison  of  Figs.  7  and  8  shows  that 
these  two  cases  are  the  same  except  that  the  positions  of  the  source  points  and  field  points  have  been 
interchanged.  In  Subsection  1.3  it  is  shown  that  if  the  detector  array  is  completely  contained  within  a 
circle  of  radius  b'  from  the  origin  and  that  every  point  in  the  source  at  a  distance  R  from  the  origin 
is  outside  of  the  Rayleigh  range  of  the  detector  array,  i.e., 

R  >>kb'2,  (36) 

then  from  Eq.  (12)  the  field  measured  by  the  detector  at  x',  can  be  given  as  the  two-dimensional 
Fourier  transform  over  the  source  distribution  by  the  Fraunhofer  approximation  of  the  third  kind, 
viz., 


^°>(x')  =  -ik/2*  ff  WRs)exp{ikR)/R]<txp{-iks  - x’)R\dsxdsy/Sl),  (37) 

s’+i^sl 


where  s  =  (sx ,  sy,  sz)  is  the  unit  vector  from  the  origin  in  the  direction  of  the  source  point  at  x. 
Upon  substitution  from  Eq.  (37)  into  Eq.  (21b)  and  by  using  Eq.  (21a)  we  find  that 

(*!’  *2)  =  (*/2t)2  ff  ff  [Hy/f.s,,  Rfr)  (exp  [ik(Rl  -  R2)]/RrR2} 
o,  d2 

x  exp  [-i*(s,  •  x[  -  Sj  •  X*)]  (Rfdsudsly/su)  (Rfa^ds^/s^),  (38) 

where  =  (su,  st),,  ru)  and  Sj  =  ($2*.  s2y,  s^)  are  the  directions  from  the  origin  to  the  source 
points  at  Xj  and  x2,  &  1  =  I  *1 1 »  &2  =  I  *2 1  >  and  D 1  and  D2  are  the  domains  of  integration  where 
su  +  Siy  ^  1  and  slx+s2y  <,  1,  respectively.  Equation  (38)  gives  the  field  cross-spectral  density 
function  between  points  xj  and  x2  within  the  detector  array  in  the  z  -  0  plane  as  a  double  two- 
dimensional  Fourier  transform  over  the  field  cross-spectral  density  function  between  the  source  posi¬ 
tions  at  R  ]S]  and  R2^2  ■  This  *s  die  second-order  form  of  the  Fraunhofer  approximation  of  the  third 
kind. 
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Fig.  8  —  Configuration  for  the  Fraunhofer 
approximation  of  the  third  kind 

We  will  consider  the  special  case  for  which  the  source  distribution  is  quasi-homogeneous  so  that 
»V(*i,  *2)  =  f[(*i  +  x2)/2] #*(*!  ~  *2).  (39> 


where  /  is  again  a  slow  function  relative  to  pi  and  n  is  a  correlation  coefficient.  Then  upon  substitu¬ 
tion  from  Eq.  (39)  into  Eq.  (38)  we  have 

«V*i-  *2)  =  (*/2t)2G(xJ  -  xi)£[(xj  +  x*)/2],  (40) 


as  shown  in  Appendix  B,  where 


and 


G(f)  =  JJ  I(Rs)  exp  (-iks  ■  ~i)R\dsxdsy/s2) 

2r 


M) 


00 


jj  n(Rs~)  exp  (- iks _  •  J)ds  _xds  _y. 


(41) 


(42) 


Reference  18,  Eq.  (2.11)  contains  a  different  derivation  of  Eq.  (40)  for  the  special  case  where  the 
source  is  completely  incoherent. 

Equation  (40)  is  different  from  Eqs.  (26)  and  (33)  in  that  it  gives  the  cross-spectral  density  func¬ 
tion  for  the  field  over  a  detector  array  if  the  sources  are  outside  of  the  Rayleigh  range  of  the  detector 
array.  The  results  in  the  earlier  sections  give  the  cross-spectral  density  function  for  detectors  outside 
of  the  Rayleigh  range  of  the  sources,  i.e.,  in  the  farfield.  These  are  very  important  differences. 
Usually  if  one  condition  applies  the  other  does  not. 

There  are  two  formal  differences  between  Eqs.  (40)  and  (33).  The  first,  and  most  important, 
difference  is  in  the  coordinates.  In  Eq.  (33)  the  Fourier  transformation  is  from  linear  coordinates 
within  the  source  plane  at  z  =  0  to  angular  coordinates  from  this  origin  to  the  field  points  far  away. 
In  Eq.  (40)  the  Fourier  transform  is  from  angular  coordinates  that  the  source  points  subtend  from  an 
origin  within  the  detector  array  within  the  z  =  0  plane  to  linear  coordinates  within  the  plane  of  the 
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detector  array.  This  is  a  very  important  difference  as  we  shall  see  when  dealing  with  inverse  prob¬ 
lems  in  the  next  subsection. 

The  second  formal  difference  between  Eqs.  (40)  and  (33)  is  the  additional  factors  appearing  in 
the  kernel  of  Eq.  (41).  In  general  R  and  sz  are  functions  of  s  and  cannot  be  removed  from  the 
integral.  However,  for  the  special  case  where  the  sources  are  all  located  on  a  hemisphere  of  a  con¬ 
stant  radius  R  centered  at  the  origin  and  the  sources  are  localized  within  a  region  that  subtends  such  a 
small  angle  from  the  origin  that  sz  can  be  approximated  by  unity,  then  Eq.  (40)  becomes 

W^\\,  x£)  =  (kR  /2t)2/ (x[  -  x£)£[(xj  -I-  x£)/2)].  (43) 

Equation  (43)  differs  from  Eq.  (33)  because  of  the  different  conditions  required  for  their  validity,  the 
location  of  the  origin  of  the  coordinates,  and  the  different  use  of  linear  and  angular  coordinates. 

2.5  Applications  of  the  Fraunhofer  Approximations  to  Inverse  Problems 

An  important  example  of  the  application  of  the  Fraunhofer  approximations  is  in  their  use  in 
inverse  problems.  In  an  inverse  problem  source  parameters  are  calculated  from  measurements  made 
of  the  field  that  they  radiate  or  scatter.  Fraunhofer  approximations  are  especially  useful  in  this  work 
since  the  relationship  between  source  and  field  parameters  appears  as  easily  inverted  Fourier 
transforms.  Leaving  the  more  complicated  case  of  the  Fraunhofer  approximation  of  the  first  kind  to 
another  study,  the  intensity  distributions  describing  the  sources  for  the  Fraunhofer  approximation  of 
the  second  and  third  kinds,  respectively,  are  given  by  inverting  Eqs.  (33)  and  (40),  as  shown  in 
Appendix  C,  to  get  the  source  intensity 

00 

Z(x')  =  /0/X2  Jj  g;(£,i»)  exp  (i •  x[)dSdii,  (44) 


if  Eq.  (29)  is  valid,  or 


00 

I(Rs)  =  G0  (cos2  6/\2R2)  Jj  g^(xL)  exp  (ifcs  ■  xL)d2x'„, 

-oo 


(45) 


if  Eq.  (36)  is  valid,  where  /0  and  G0  are  constants, 

**(*!  -  *2.  y[  -  yi)  =  ^(*1.  *$)/>/ i*V(x{,  *2>. 

iffi  x'\  +  y'\  <  b’2  and  x'\  +  y'\  s  b’2, 

~  0  otherwise  (46) 

is  the  normalized  correlation  coefficient  measured  with  the  detector  array  over  the  z  =  0  plane,  and 

lor  ~s2x’sl y  ~  $2y)  —  R,^(H]sl.f?2S2)exP l— *^(^1  —  ^lsl)  ^V(^2*2» ^2®2)> 

iffi  sfx  +  s2y  rs  1  and  six  +  sly  S  1, 

4  0  otherwise,  (47) 
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is  a  normalized,  phase-shifted  correlation  coefficient  measured  over  the  detector  array  away  from 
z  =  0  plane.  In  Eqs.  (46)  and  (47)  the  correlation  coefficients  are  spatially  stationary. 

Note  that  the  main  formal  difference  between  Eqs.  (44)  and  (45)  is  the  different  coordinates  in 
the  two  Fourier  conjugate  spaces  and  the  obliquity  factors.  In  Eq.  (44)  the  transform  is  from  a  space 
linear  with  the  difference  in  angle  subtended  from  the  origin  at  the  sources  by  the  points  at  which  the 
field  correlation  function  is  measured  into  a  space  linear  with  the  source  plane  coordinates,  as  shown 
in  Fig.  7.  In  Eq.  (45)  the  transform  is  from  coordinates  linear  with  the  position  of  the  detectors  in 
the  z  =  0  plane  into  angular  coordinates  linear  with  the  angle  that  the  source  point  subtends  from  the 
+z  axis.  Since  the  transformation  from  linear  to  angular  coordinates  is  not  a  linear  transformation 
(see  Eqs.  (15)  and  (16)),  this  difference  is  important  in  the  processing  of  field  correlation  function 
data  to  obtain  source  intensity. 

An  important  example  of  this  inverse  problem  is  radio  astronomy  where  correlations  of  the  field 
radiated  by  astronomical  sources  are  measured  by  an  array  of  antennas  at  the  surface  of  Earth  and 
then  are  used  to  calculate  the  secondary  source  distribution  over  a  hemisphere  outside  of  the  Rayleigh 
range  of  the  array.  Equation  (45)  is  the  basis  for  the  present  radio  telescopes  that  image  sources  out¬ 
side  of  the  Rayleigh  range  of  the  antenna  array  [18].  Equation  (44)  could  be  used  as  the  basis  of  a 
new  application  of  a  radio  telescope  to  image  a  compact  source  distribution  that  is  within  the  Rayleigh 
range  of  the  antenna  array  but  for  which  the  antenna  array  is  in  the  farfield  of  the  source.  However, 
because  of  the  formal  differences  between  Eqs.  (44)  and  (45)  the  data  from  a  radio  telescope,  used  to 
image  a  small  source  distribution  inside  the  Rayleigh  range  of  the  antenna  array,  will  require  different 
processing  to  image  sources  outside  of  the  Rayleigh  range  of  the  array.  Because  the  coordinate 
transformation  from  linear  to  angular  coordinates,  required  in  Eq.  (44)  before  Fourier  transformation, 
is  dependent  on  R(x+)  (the  range  from  each  antenna  to  the  origin  within  the  source  distribution),  the 
resulting  image  will  be  accurate  only  for  the  sources  that  are  near  enough  to  the  origin  (in  agreement 
with  the  condition  in  Eq.  (29)).  It  therefore  appears  that  this  transformation  of  the  correlation  func¬ 
tion  data  to  angular  coordinates  before  the  Fourier  transformation  refocuses  the  telescope  on  the  local¬ 
ized  source  distribution  and  is  analogous  to  focusing  a  camera  on  an  object  that  is  too  close  to  be  in 
focus  when  the  camera  is  focused  at  infinity.  This  procedure  is  discussed  in  Section  3  (see  also  Ref. 
19). 


2.6  Differences  Between  Fraunhofer  Approximations 

This  subsection  summarizes  and  discusses  some  of  the  more  important  differences  discovered 
between  the  Fraunhofer  approximations  of  the  first,  second,  and  third  kinds. 

A  major  mathematical  difference  between  the  Fraunhofer  approximation  of  the  first  kind  (Eq. 
(24))  and  the  others  is  that  it  represents  a  Fourier  transform  from  a  six-dimensional  space  to  a  four¬ 
dimensional  space.  This  is  because  xj  and  x£  in  Eq.  (24)  are  ordinary  vectors  with  three  independent 
components  each,  while  s,  and  Sj  are  unit  vectors  with  only  two  independent  components  each.  It  is 
therefore  not  clear  that  Eq.  (24)  can  be  generally  inverted  to  apply  to  inverse  problems.  It  appears 
likely,  however,  that  inversion  is  possible  for  some  special  cases.  The  Fraunhofer  approximations  of 
the  second  and  third  kind  in  Eqs.  (31)  and  (38)  are  not  as  difficult  to  invert.  Appendix  C  shows  how 
this  can  be  done  for  quasi-homogeneous  sources  when  certain  approximations  are  made. 

The  important  differences  in  the  Fraunhofer  approximations  of  the  second  and  third  kinds  are 
both  in  form  and  in  conditions  for  validity.  The  difference  in  conditions  for  validity  were  described 
during  derivation.  The  Fraunhofer  approximation  of  the  second  kind  is  valid  only  if  the  domain  of 
the  source  distribution  is  small  enough  that  the  detectors  can  all  be  located  outside  of  the  Rayleigh 
range  of  the  source  distribution.  On  the  other  hand,  the  Fraunhofer  approximation  of  the  third  kind  is 
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valid  only  if  the  detector  array  is  small  enough  that  the  sources  can  all  be  placed  outside  of  the  Ray¬ 
leigh  range  of  the  detector  array.  Usually  if  one  approximation  aoplies  the  other  does  not.  In  cases 
where  both  approximations  apply  the  formal  differences  between  them  can  be  removed  by  small  angle 
approximations  (see  Subsection  1.4). 

The  formal  differences  between  the  Fraunhofer  approximations  of  the  second  and  third  kind  are 
twofold,  differences  in  both  the  coordinate  systems  and  the  obliquity  factors.  It  has  been  pointed  out 
(Eq.  (16))  that  the  coordinate  transformations  from  linear  to  angular  coordinates  in  Eqs.  (38),  (41), 
(43),  or  (44)  needed  to  prepare  the  data  for  Fourier  transformation  are  in  general  nonlinear  transfor¬ 
mations  and  do  not  commute  with  the  Fourier  transformation.  No  such  coordinate  transformations 
are  required  for  Eqs.  (24),  (26),  (31),  (33),  or  (45),  for  which  the  transformed  data  come  out 
automatically  in  angular  coordinates  and  need  to  be  transformed  into  linear  coordinates  only  if 
desired.  Thus  the  formal  differences  in  processing  data  by  using  different  Fraunhofer  approximations 
is  nontrivial  (viz.  Section  3). 

2.7  Conclusions 

The  differences  between  Fraunhofer  approximations  of  the  first,  second,  and  third  kinds,  which 
are  discussed  in  terms  of  scalar  diffraction  theory  in  Section  1,  are  also  found  to  exist  when  these 
theories  are  extended  to  describe  the  propagation  of  the  cross-spectral  density  function.  The 
Fraunhofer  approximation  of  the  third  kind  leads  to  the  formulas  that  are  usually  used  in  radio  astron¬ 
omy  to  process  the  data  from  a  radio  telescope  to  image  sources  in  the  farfield  of  the  antenna  array. 
The  use  of  the  Fraunhofer  approximation  of  the  second  kind,  instead,  leads  to  modified  formulas.  It 
appears  that  these  formulas  show  how  to  process  these  data  differently  so  that  the  radio  telescope  is 
refocused  to  image  objects  in  the  nearfield  of  the  antenna  array. 

3  REFOCUSING  A  RADIO  TELESCOPE  TO  IMAGE  SOURCES  IN  THE 
NEARFIELD  OF  THE  ANTENNA  ARRAY 

3.1  Introduction 

Radio  telescopes  have  been  used  since  the  1950s  to  image  radio  sources  that  are  very  far  from 
the  antenna  arrays.  There  has  been  some  interest  at  the  Naval  Research  Laboratory  in  building  new 
radio  telescopes  with  very  large  space  borne  antenna  arrays  to  improve  resolution.  Some  of  these 
arrays  would  be  so  large  that  some  sources  within  the  solar  system  would  be  within  their  nearfield  so 
that  the  usual  methods  for  forming  images  fail.  Some  experimental  work  has  already  been  done 
toward  refocusing  a  radio  telescope  to  look  at  nearby  objects  but  without  the  support  of  any  detailed 
theoretical  understanding.  This  work  has  been  done  principally  by  use  of  optical  analogies. 

In  the  present  section  a  theoretical  foundation  is  provided  for  understanding  the  refocusing  prob¬ 
lem  in  a  more  systematic  way  by  use  of  optical  coherence  theory.  The  first  subsection  describes  the 
operation  of  a  conventional  radio  telescope  focused  at  infinity.  The  second  subsection  extends  this 
theory  to  describe  some  methods  for  refocusing  the  telescope  to  image  objects  that  are  much  closer  to 
the  antenna  array.  This  study  is  based  on  the  formulation  derived  in  Sections  1  and  2  (also  see  Refs. 
15  and  20). 

3.2  Conventional  Farfield  Imaging 

Conventionally  a  radio  telescope  is  used  to  image  radio  sources  that  are  in  the  farfield  of  the 
antenna  array  of  the  instrument.  This  is  analogous  to  a  camera  focused  at  infinity.  If  the  radiation 
source  is  in  the  farfield  of  the  antenna  array,  it  makes  no  difference  to  the  adjustment  of  the  telescope 
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or  to  the  computation  of  the  image  just  how  far  the  source  is  from  the  telescope.  All  objects  are  in 
focus  simultaneously  if  they  are  in  the  farfield  of  the  antenna  array.  It  is  important  to  remember  tha» 
in  radio  astronomy  the  term  farfield  refers  to  the  source  being  in  the  farfield  of  the  antenna  array.  As 
we  shall  soon  see  this  is  quite  different  than  having  the  antenna  in  the  farfield  of  the  sources  just  as 
discussed  in  Sections  1  and  2.  The  conventional  method  by  which  a  radio  telescope  images  radio 
sources  in  its  farfield  is  described  in  this  subsection  for  comparison  to  the  new  technique  for  refocus¬ 
ing  the  telescope  to  image  sources  in  its  nearfield. 

Consider  an  antenna  array  for  a  radio  telescope  arranged  over  the  z  =  0  plane  near  the  origin, 
as  shown  in  Fig.  9,  irradiated  by  radio  waves  from  sources  in  its  farfield.  Since  the  sources  are  all  in 
the  farfield  of  the  antenna  array,  to  image  them  we  need  only  calculate  the  intensity  over  a  hemi¬ 
sphere  of  radius  R  from  the  origin. 


Fig.  9  —  Conventional  radio  telescope  focused  at  infinity 

The  antennas  measure  the  amplitude  of  one  component  of  the  transverse  electric  field  at  points 
indicated  by  the  position  vector  x'  from  the  origin  as  shown  in  Fig.  9.  The  signal  from  each  antenna 
is  then  processed  as  shown  in  Fig.  10.  The  signals  are  detected  by  a  receiver,  reduced  to  a  con¬ 
venient  IF,  and  bandlimited  to  a  narrow  frequency  spectrum.  Then  the  signals  are  passed  through 
delay  lines  that  point  the  array  toward  the  source  so  that  the  normal  to  the  plane  of  the  antenna  array 
points  in  the  direction  of  interest.  The  resulting  signals  are  then  multiplied  together  and  the  product 
time  averaged  to  form  correlation  functions  between  signals  from  antenna  pairs,  viz., 

r(x,,  x2,  0)  =  lim  1/27  1  *(X[,  f)  ♦*(x2,  t)dt ,  (48) 

where  ♦(x,  t)  represents  the  complex  field  amplitude  of  the  traveling  wave  propagating  from  the 
source  to  the  antenna  array,  after  it  has  been  detected  and  bandlimited  by  the  telescope. 

We  represent  the  second-order  statistical  properties  of  the  field  by  use  of  either  the  mutual 
coherence  function  or  the  cross-spectral  density  function  of  optical  coherence  theory  (for  general 
description  see  Refs.  3  and  21,  and  for  applications  to  radio  astronomy  see  Ref.  18).  The  second- 
order  statistical  properties  of  the  field  are  completely  specified  by  the  mutual  coherence  function 

T(x,,  x2,  r)  =  lim  1/27  f  *(x,,  t  +  t)  ¥’(x2,  t)dt ,  (49) 

r— od  J -t 
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which  is  usually  defined  by  using  an  ensemble  average  instead  of  the  time  average  used  in  radio 
astronomy.  These  are  equivalent  if  the  statistics  of  the  source  are  assumed  to  be  ergodic.  Note  that 
Eq.  (49)  is  a  simple  extension  of  the  function  described  by  Eq.  (48),  which  is  measured  by  the  radio 
telescope.  The  cross-spectral  density  function  is  then  defined  by  the  equation 

Wu(xu  x2)  =  |°°  r(x„  x2,  r)  e'“T  dr,  (50) 

as  the  Fourier  transform  of  the  mutual  coherence  function. 

The  procedure  for  forming  an  image  depends  critically  on  some  important  assumptions.  First 
we  must  assume  that  the  antenna  array  is  bounded  so  that  all  antennas  are  within  a  circle  of  radius  a 
from  the  origin  within  the  z  =  0  plane.  Second,  we  must  assume  that  all  of  the  sources  are  in  the 
farfield  of  the  antenna  array  so  that 

R  >>  ka2,  (51) 

where  k  =  oj/c  .  And  third,  we  assume  that  the  sources  of  the  field  are  quasi-homogeneous  [17]. 
Then  it  is  Well  known  (see  Eq.  (45)  and  the  preceding  derivation)  that  the  intensity  associated  with  the 
measured  component  of  the  electric  field  over  the  hemisphere  at  infinity  (or  equivalently  at  fixed  dis¬ 
tance  R  from  the  origin)  is  given  by  the  Fourier  transform  of  the  normalized  cross-spectral  density 
function,  i.e., 

00 

I(Rs)  =  G0  cos2  0/X2/?2  ||  g(u)e~ika,d2u,  (52) 
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where 


S(xj  -  x£)  4  WJx{,  xj)^,  x£), 

iffi  j X|  |  s  a  and  | x2  |  s  a, 

=  0  otherwise,  (53) 


is  a  spatially  stationary  correlation  coefficient,  s  is  the  unit  vector  pointing  from  the  origin  in  the 
direction  of  an  intensity  point  on  the  hemisphere,  6  is  the  angle  that  s  makes  with  the  +xaxis, 
u  =  («,  v)  where 


u  =  x{  - 
and 

v  =  y[  -  y'z .  (54) 

are  difference  coordinates  between  antenna  pairs  on  the  z  =  0  plane,  and  G0  is  a  constant.  If  some 
important  restrictions  are  made  on  the  operation  of  the  radio  telescope  (see  Ref.  18,  Section  II  for  a 
detailed  discussion),  we  can  approximate  the  cross-spectral  density  function  in  Eq.  (52)  by  the  mea¬ 
sured  correlation  function  defined  by  Eq.  (48),  i.e.,  we  can  use 

g(u)  =  r(x|,  x£,  0)/  Vr(x|,  xj,  0)  T(x£,  x£,  0) 
iffi  |  x(  |  s  a  and  |  x£ |  &  a, 

—  0  otherwise,  (55) 

in  the  kernel  of  Eq.  (52),  as  shown  in  Appendix  D.  Therefore,  the  image  may  be  obtained  by  using 
Eq.  (52)  by  normalizing  the  data  and  then  performing  a  Fourier  transform  from  the  two-dimensional 
«v  plane  of  antenna  position  differences  to  the  two-dimensional  angular  space  given  by 

sx  =  x/R 


and 


sy  =  y/R,  (56) 

where  x  is  a  position  vector  from  the  origin  to  a  point  on  the  hemisphere  as  shown  in  Fig.  9. 

If  the  sources  are  nearer  to  the  antennas  so  that  they  are  not  in  the  farfield  of  the  array,  then, 
under  some  circumstances,  there  is  still  a  way  to  refocus  the  telescope  to  image  these  sources. 

3.3  NearfWd  Imaging 

Consider  Fig.  11,  which  shows  a  localized,  planar,  secondary,  quasi-homogeneous  source  distri¬ 
bution  over  the  z  =  0  plane  radiating  radio  waves  toward  an  antenna  array  of  a  radio  telescope  that  is 
located  in  the  z  =  D  plane. 
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We  assume  that  the  sources  are  not  in  the  farfield  of  the  antenna  array  so  the  D  is  too  small  to 
satisfy  Eq.  (51).  Even  so  a  method  exists  to  form  an  image  by  using  the  same  radio  telescope 
described  schematically  in  Fig.  10  but  slightly  modifying  the  processing  procedure  to  refocus  the  tele¬ 
scope  onto  the  nearby  sources. 

This  process  requires  some  important  assumptions.  First,  we  assume  that  the  sources  are 
bounded  so  that  they  are  all  contained  within  a  circle  of  radius  b  from  the  origin  in  the  z  =0  plane. 
Second  we  assume  that  the  antenna  array  is  in  the  farfield  of  the  sources  (it  is  the  other  way  around 
for  conventional  farfield  imaging)  so  that 


R  >>  kb2. 


Then  as  shown  in  Eq.  (44),  the  intensity  over  the  z  =  0  plane  is  given  by 


/(*')  A  /0/X2  |J  g  ’(£)  exp  (ikx'  ■  &d% 


(57) 


(58) 


where 


g'(x,/R,  -  x2/R2)  A  ww(x ;,  x£)  exp  [~ik(R ,  -  R2 ]/  y/WJx[,  x[)Wu(xi,  x$), 
iffi  |  X,  |  /J?  s  1  and  |*jl/R  £  1, 

A  0  otherwise,  (59) 

must  be  calculated  from  the  experimental  correlation  data. 
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Inspection  of  Eqs.  (58)  and  (59)  suggests  a  straightforward  computational  method  for  refocusing 
the  telescope  to  image  nearby  objects.  As  shown  in  Appendix  E,  under  almost  the  same  conditions 
used  with  the  telescope  focused  at  infinity,  the  correlation  coefficient  in  Eq.  (59)  can  be  approximated 
by 


g'(x,//?,,  x2/R2)  =  r(x„  x2,  0)  exp  [~ik(Ri  -  /?2)]/Vr(x1,  x,.  0)r(x2,  x2,  0), 
iffi  xf  +  y2  <  /?2  and*2  +  y2  <  R2  , 

=  0  otherwise,  (60) 

by  using  the  measured  correlation  function  in  Eq.  (48).  Then  to  form  an  image  of  the  fields  over  the 
hemisphere,  g'  must  be  transformed  from  linear  u  ,v  coordinates  within  the  antenna  plane  to  the  new, 
angular  coordinates  given  by 

«  =  x,/Vx?  +  y\  +  D2  -  x2/Vx|  +  y f  +  D2, 

V  =  yi/Vx?  +  y?  +  D2  -  y|/Vx|  +  y\  +  D2,  (61) 

and  then  Fourier  transformed.  Equation  (61)  is  a  nonlinear  coordinate  transformation  that  has  a  non¬ 
trivial  effect  on  the  data  after  Fourier  transformation. 

Inspection  of  Eqs.  (58)  and  (59)  suggests  that  it  might  perhaps  be  useful  to  refocus  the  telescope 
electronically.  This  could  be  done  by  adjusting  the  delay  lines  in  Fig.  10  so  that  the  antenna  array  is 

electronically  moved  from  the  plane  to  a  hemisphere  of  radius  R  from  the  origin  as  shown  in  Fig.  1 1 . 

This  can  be  done  since  the  antennas  are  in  the  farfield  of  the  sources  and  the  farfield  amplitude  is 
constant  along  any  direction  from  the  origin  except  for  at  IR  decay  and  the  usual  exp  (ikR)  phase 
shift.  This  use  of  the  delay  lines  has  the  effect  of  refocusing  the  antenna  array  onto  the  sources  in 
analogy  to  a  refocusing  a  camera  on  a  nearby  object.  With  refocusing  the  function  in  Eq.  (59) 
becomes  the  same  normalized  correlation  coefficient  given  in  Eq.  (53)  and  it  can  be  calculated  from 
the  measured  correlation  function  by  using  Eq.  (55)  just  as  for  the  case  when  the  telescope  was 
focused  at  infinity. 

If  we  adjust  the  delay  lines  as  suggested  so  that  in  Eq.  (58)  Rx  =  R2  =  D  a  constant,  then  Eq. 
(61)  reduces  to 


$  =  u/D 


and 


n  =  v/d, 


(62) 


so  that  Eq.  (58)  becomes 


/(*')  =  /0A2D2  JJ  g(u)  exp  (ilex' •  a /D2)dudv,  (63) 
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which  is  essentially  the  same  as  Eq.  (52)  except  for  the  trivial  cos  6  factor  and  the  fact  that  the  u  ,v 
coordinates  in  Eq.  (63)  correspond  to  the  antenna  position  differences  after  they  are  projected  onto  the 
hemisphere,  not  those  on  the  actual  physical  antenna  plane.  As  shown  in  Appendix  E,  the  correlation 
coefficient  in  the  kernel  can  be  approximated  by  using  the  measured  correlation  function  as  given  by 
Eq.  (60)  but  with  exp  (—ik(Rl  -  R2))  =  1.  Thus  the  data  can  be  processed  essentially  as  before 
when  the  telescope  was  focused  at  infinity. 

In  appears  clear  from  these  results  that  we  have  to  either  refocus  the  telescope  physically  by  use 
of  the  phase  shifters  or  refocus  it  mathematically  by  use  of  a  modified  data-processing  procedure.  In 
either  case  we  should  be  able  to  form  an  image  of  the  nearby  object  in  a  similar  manner  to  an  object 
in  the  farfield. 

3.4  Conclusions 

It  appears  possible  to  refocus  a  radio  telescope  to  image  radio  sources  within  the  nearfield  of  the 
antenna  array  if  the  antenna  array  is  in  the  farfield  of  the  source  distribution.  This  is  done  by 
refocusing  the  telescope  onto  the  sources  in  analogy  to  a  photographic  camera.  Refocusing  can  be 
done  either  electronically  by  curving  the  antenna  array  into  a  portion  of  a  hemisphere  with  a  center  at 
the  source,  or  by  processing  the  correlation  function  data  obtained  from  the  antennas  in  a  somewhat 
modified  manner.  Neither  of  these  techniques  would  require  any  modification  to  most  existing  tele¬ 
scopes.  The  modifications  to  the  data-processing  software  should  not  be  extensive. 
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Appendix  A 

RECIPROCITY  THEOREM 


The  reciprocity  theorem  used  in  antenna  theory  treats  antennas  as  circuit  elements  and  usually 
relates  only  the  input  and  output  currents  and  voltages  for  coupled  antennas.  Some  work  has  been 
done  to  obtain  reciprocity  theorems  for  electromagnetic  fields  [Al,  Chapter  9]  but  this  work  is  diffi¬ 
cult  to  relate  to  field  propagation  as  discussed  here.  Therefore,  a  generalized  field  reciprocity 
theorem  is  derived  in  this  appendix  in  the  form  needed  to  obtain  Eq.  (12)  from  Eq.  (10). 

First,  consider  the  outgoing  field  in  Fig.  2.  Assume  that  the  field  arises  from  primary  sources 
in  the  z  <  0  half  space  and  propagates  to  the  right  into  the  z  2s  0  half  space.  Let  the  z  2s  0  half 
space  contain  only  a  homogeneous,  isotropic,  source-free  medium.  We  also  assume  that  the  sources 
are  far  enough  from  the  z  =  0  plane  that  the  effects  of  evanescent  plane  waves  on  the  fields  in  the 
z  >  0  half  space  can  be  neglected.  Then  we  can  always  expand  the  scalar  amplitude  of  the  field  in 
this  half  space  into  an  angular  spectrum  of  homogeneous  plane  waves  in  the  form 


W*2)  =  JJ  A,(/>2.  92)  exp  [i*(/>2*2  +  92),2  +  m222)]4>2^92. 


pi  +«2  s  1 


where 


m2  =  Vl 


-Pi  ~ 


Pi  +  ql  ^  1. 


(Al) 

(A2) 


The  angular  spectrum  is  found  by  taking  the  Fourier  transform  of  Eq.  (Al)  after  setting  z  =  0  to  get 


A i(P2>  92)  =  1A2  H  ^](0)(x2)  exp  [-!*(/> 2*2  +  92y2)]d!r2Ay2,  (A3) 

—  00 


and  the  field  over  the  hemisphere  of  radius  R  from  the  origin  is  found  by  integrating  Eq.  (Al)  by  use 
of  the  method  of  stationary  phase  to  get 

4>\a‘\^i) - i\z2/R  A\(x2/R,  y2/R)eikR  /R,  (A4) 

asymptotically  as  kR  becomes  large. 

Next  consider  the  reciprocal  field  illustrated  in  Fig.  3.  Assume  that  this  field  arises  from  pri¬ 
mary  sources  to  the  right  of  the  hemisphere  in  Fig.  3  and  propagates  to  the  left  toward  the  z  =  0 
plane.  We  assume  that  the  region  to  the  left  of  the  hemisphere  contains  only  a  homogeneous,  isotro¬ 
pic,  source-free  medium.  We  also  assume  that  the  primary  sources  are  far  enough  away  from  the 
hemisphere  that  we  can  neglect  the  effects  of  evanescent  plane  waves  on  the  field  to  the  left  of  the 
hemisphere.  Then  over  the  region  between  the  z  =  0  plane  and  the  hemisphere  we  can  always 
expand  the  scalar  amplitude  of  the  field  into  an  angular  spectrum  of  homogeneous  plane  waves  in  the 
form 

=  JJ  A2(/?3,  q3)  exp  [ik(pyx3  +  y3  -  m 323)^3^ 3,  (A5) 

Pi  +?3  S’ 
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where 

m3  =  Vl  -  pi  +  q\  ,  p\  +  q\  S  1.  (A6) 

The  angular  spectrum  is  obtained  from  Eq.  (A5)  as 

GO 

^2(P3*  $3 )  =  1  A2  Jf  ^j0)(* 3)  exp  [~ik(psx $  +  q&QWtfyi  (A7) 

—  OO 

and  the  secondary  source  distribution  over  the  hemisphere  is  found  from  Eq.  (A5)  to  be 

~  i\z2/RA2(-x3/R,  -  yt/Rye-^/R,  (A8) 

asymptotically  as  kR  becomes  large. 

Let  the  outgoing  field  in  Fig.  2  be  called  field  one  and  the  incoming  field  in  Fig.  3  be  called 
field  two.  We  assume  that  field  one  has  an  arbitrary  secondary  source  distribution  given  by 

*j0)(x$)  £  G(xi),  (A9) 

over  the  z  =  0  plane.  Then  by  substitution  from  Eq.  (A9)  into  Eq.  (A3)  we  obtain  the  angular  spec¬ 
trum 

00 

A i(/*2>  92)  =  G(p2,  q2 )  =  1/A2  Jf  G(xi)  exp  [-ik(pjx2  +  92? 2)!^ 2*^2-  (A10) 

-OO 

By  substitution  from  Eq.  (A10)  into  Eq.  (A4)  we  find  the  field  distribution  over  the  hemisphere  to  be 

~  -i\z2/RG(x2/R,y2/R)eikR/R.  (All) 

We  now  define  field  two  in  Fig.  3  to  be  the  reciprocal  of  field  one  by  requiring  that  the  secondary 
source  distribution  over  the  hemisphere  in  Fig.  3  be  the  complex  conjugate  of  the  field  distribution 
given  in  Eq.  (All)  for  field  one  over  the  identical  hemisphere  in  Fig.  2,  i.e., 

^O0)(*3)  ^  =  i\zi/RG*(xi/R,y2/R)e-ikR/R.  (A12) 

By  comparing  Eq.  (A  12)  with  Eq.  (A8)  it  is  clear  that  field  two  must  have  the  angular  spectrum 

A2(p2,  93)  =  G*(—p2,  -qj).  (A  13) 

Upon  substituting  from  Eq.  (A10)  into  Eq.  (Al)  we  find  that 

^i(*2>  =  G(P2>  92)  exp  [ik(p2X2  +  q&2  +  m^dpjdq^  (A14) 

Pi  +?|s> 
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and  then  by  substituting  from  Eq.  (A  13)  into  Eq.  (AS),  taking  the  complex  conjugate,  setting 
Pi  =  -pj  and  q3  =  -q3,  we  find  that 

^2*(*3)  =  fj  G(Pi,  q3)  exp  [ik(piX3  +  +  m^i))dpydqy  (A15) 

Pi  S 1 

By  comparing  Eqs.  (A14)  and  (A15)  it  is  obvious  that  at  equivalent  points  in  the  two  spaces, 
where  x2  =  x3  =  x,  the  fields  take  on  values  that  are  complex  conjugate  to  each  other  so  that 

iM*)  =  &*(x).  (A16) 

Thus  if  we  derive  a  forward  propagation  integral  for  field  one  in  the  form 

OB 

*,(*")  =  Jf  *V2x\  (A  17) 


from  any  plane  z  =  z'  to  any  later  field  point  x",  both  contained  within  the  region  between  the 
z  =  0  plane  and  hemisphere  in  Fig.  2,  then  by  substitution  from  Eq.  (A16)  into  Eq.  (A17)  and  taking 
the  complex  conjugate  we  have 


hi*")  *  If  iM*')*V',x>/2x\  (A  18) 


for  the  inverse  propagator  for  field  two  from  the  corresponding  z  =  z'  plane  back  to  the  earlier  field 
point  at  x". 

Thus  we  have  proved  a  useful  reciprocity  theorem  that  can  be  summarized  as  follows:  If  two 
electromagnetic  fields,  propagating  in  opposite  directions,  taken  one  at  a  time,  within  a  source-free, 
isotropic,  and  homogeneous  medium,  are  composed  only  of  homogeneous  plane  waves  such  that  their 
angular  spectra  are  related  by 


A,(p,  q)  £  Al(-p,  -q), 

then  they  are  reciprocal  so  that  the  forward  propagator  for  one  field  (i.e.,  the  kernel  h  in  Eq.  (A17)) 
is  the  complex  conjugate  of  the  inverse  propagator  (i.e.,  the  kernel  h*  in  Eq.  (A  1 8))  for  the  other 
field. 
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Appendix  B 

DERIVATION  OF  EQUATION  (40) 


Substituting  from  Eq.  (39)  into  Eq.  (38)  and  transforming  variables  by  using 


S_  —  S(  —  S2, 

S+  =  (S!  +  S2)/2, 

=  S+  +  s_/2. 

S2  =  s+  —  s_/2, 

(Bl) 

xl  =  xj  -  x2, 

x'+  =  (xj  +  xg)/2. 

xj  =  x;  +  x'_/2, 

x2  =  x\  -  xl/2, 

(B2) 

then  defining  R  +  and  R  _  such  that 

xt  =  /?i  ss  =  R+s+  +  /?_  s_/2, 

*2  =  ^2*2  *  K +s+  ~  R-S-/2,  (B3) 


we  obtain 


QO  OO 

W^x i,  xg)  =  (k/2r)2  JJ  Jf  (/(i?+s+)At(i?_s_)  [exp  [«*(*,  -  *2)]]| 

—  00  —  oo 

x  exp  [-iifc(s+  •  xl  +  s_  •  x+)](/? lR2ds  +xds  +y /s izS2z)ds _xds (B4) 

We  must  assume  that  /(x+)  =  0  outside  of  the  domains  where  S)  and  Sj  point  into  the  z  >  0  half 
space  so  that  we  can  extend  the  bounds  of  integration  in  Eq.  (38)  to  infinity  before  we  transform 
coordinates  to  get  Eq.  (B4). 

We  now  make  use  again  of  our  previous  assumption  that  /t(x_)  =  0  outside  of  the  domain 
where  |x_  |  s  d.  First,  we  can  assume  that  d  is  small  enough  that  /? i(xt),  R2(x 2),  stz  and  s^,  like 
7(x+),  do  not  vary  appreciably  over  the  domains  where  X|  and  x2  vary  by  d,  so  that  we  have  the  rela¬ 
tions 


/?,=/?+  A  R, 
R2~  R+  A  R, 
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and 


SU  ~  s2z  -  sz>  (B5) 

where  the  geometrical  significance  of  these  parameters  are  shown  in  Fig.  Bl.  Second,  we  assume 
that  because  of  the  small  d,  the  angle  between  st  and  $2  in  Fig.  Bl  is  small  enough,  for  all  values  of 
X|  -  x2  for  which  the  kernel  of  the  integral  in  Eq.  (B4)  does  not  vanish,  that  we  can  treat  s+  as  a 
unit  vector  s. 


By  substituting  from  Eq.  (B5)  into  Eq.  (B4)  and  separating  the  independent  integrations  we  find 
that 


00 

H^(xJ,  x2)  =  (k /2x)2  JJ  I(Rs)  exp  (-iks  ■  xL )R2(dsxdsy/s*) 


oo 

X  /i(/?_s_)  exp  (-jfcs_  •  x'+)ds -xds -y. 
—  00 


(B6) 


Now,  since  I(Rs)  =  0  outside  of  the  domain  where  (sf  +  sf)  <  1,  we  can  reset  the  range  of 
integration  over  s  so  that  Eq.  (B6)  becomes  the  same  as  Eq.  (40)  with  the  two-dimensional  Fourier 
transforms  defined  as  given  by  Eqs.  (41)  and  (42). 


Appendix  C 

DERIVATION  OF  EQUATIONS  (44)  AND  (45) 


To  derive  Eq.  (45)  from  Eqs.  (40)  and  (41)  we  first  normalize  the  correlation  function  in  Eq. 
(40)  to  form  the  correlation  coefficient  defined  by 

M*(x{,  *2)  =  W£\[,  xft/yJW^x [,  x[ )W+(xi,  xft,  (Cl) 

to  obtain 

Mx{,  *J)  =  G(x|  -  xi)Ml(x I  +  x^)/2]/  [g(0)  V*(x{)ji(x$)]  .  (C2) 


Since  /(x+)  in  Eq.  (39)  (as  well  as  Rlt  R2,  slz,  and  s^)  is  a  slow  function  relative  to  /t(x_),  we 
know  from  Fourier  theory  that  £(x'+)  is  a  slow  function  relative  to  G(x'_)  in  Eq.  (40).  Thus  we  can 
assume  that  G(jtj  -  x'2)  —  0  except  over  some  domain  where  |xj  -  x£|  ^  d',  and  that  within  any 
domain  where  x|  and  x2  differ  by  no  more  than  d '  we  have 

A(x i)  =  m(*2)  ~  AK*!  +  xi)/2], 
so  that,  upon  substitution  from  Eq.  (C3)  into  Eq.  (C2)  we  have 

^(x{,  x£)  k  g+(x[  -  x2)  =  G(xJ  -  x2)/G0, 

where  G0  =  G(0)  is  a  constant. 

Since  the  source  distribution  is  defined  physically  only  for  points  in  directions  from  the  origin  to 
the  z  >  0  half  space  we  can  define  I(Rs)  in  Eq.  (41)  to  have  a  value  of  zero  for  (s?  +  Sy)  >  1. 
Thus  we  can  define  the  Fourier  inverse  of  Eq.  (41)  in  the  usual  way.  Upon  substitution  from  Eq. 
(C4)  into  this  inverse,  we  obtain  Eq.  (45). 

The  derivation  of  Eq.  (44)  proceeds  in  a  similar  manner.  Upon  substitution  from  Eq.  (33)  into 
Eq.  (Cl)  we  have 

M*(*i S|,  R2^2>  =  /(*  1  -  s2)  A  [(Si  +  S2)/2]  exp  [/A (/? ,  -  R2)/  [/(0)  V^s,)^)]  .  (C5) 


(C3) 

(C4) 


Because  ft  is  a  slow  function  relative  to  /  we  have 

A(si)  ~  Ato)  ~  AKsi  +  S2)/2]-  (C6) 

within  the  domain  where  |  st  -  S2I  s  d/R,  so  that  Eq.  (C5)  becomes 

M*iS„  R2S2)  =  /( Sj  —  S2)//o  exp  [ik(R\  —  If  2)],  (C7) 

where  70  =  1(0)  is  a  constant. 
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Equation  (Cl)  has  physical  meaning  only  for  unit  vectors  s,,  and  S2,  that  point  from  the  origin 
to  the  7  sO  half  space;  therefore,  if  we  define 

-  S 2x,  Sly  -  S 2y )  =  R&)  «P  I~ *A(Jf j  -  f?2)]> 

iffi  +  jfy  <  1  and  s|,  +  s%y  s  1, 

=  0  otherwise,  (C8) 

substitute  Eq.  (C8)  into  Eq.  (Cl),  solve  for  / ,  and  then  substitute  the  result  into  the  Fourier  inverse 
of  Eq.  (34)  we  obtain  Eq.  (44). 


Appendix  D 

CALCULATION  LEADING  TO  EQUATION  (55) 


The  cross-spectral  density  function  over  the  antenna  plane  in  Fig.  9  is  given  by 


K0)(x{,  x$)  =  (Jt/2x)2G(w)(xJ  -  x$.  to)£(eB)[(xi  +  x*)/2,  to],  (Dl) 


where 


and 


G(00)(x'_,  to)  ^  J|  #2/<00>(/?s)/cos2  6  exp  (iks  •  x'_)dsxdsy,  (D2) 


fi(a,)(x+,  a)  =  jj  M(“>(/?s)  exp  (iks  •  x'+)dsxdsy  (D3) 


(see  Eqs.  (40),  (41),  and  (42)),  /<00)(x)  is  the  intensity  distribution  over  the  hemisphere  in  Fig.  9,  and 
jt<oo)(x_)  is  the  complex  degree  of  spectral  coherence  over  the  hemisphere.  The  field  over  the  hemi¬ 
sphere  is  assumed  to  be  quasi-homogeneous  so  that  /(M)(x)  is  a  slow  function  relative  to  /QO)(x_), 
i.e.,  ^<0D)(x1  -  x2)  =  0  except  within  the  domain  where  |xj  -  x2|  <  d,  and  /<00>((x,  +  x2)/2)  is 
essentially  constant  over  any  domain  in  which  |  (X|  +  x2)  |  /2  changes  by  no  more  than  d. 

Upon  substitution  from  Eq.  (Dl)  into  the  Fourier  inverse  of  Eq.  (50)  we  find  that 
r(0,(x|,  x2,  r)  =  l/(2x)  f  [(fc/2x)2d'“)(xi  -  x$,  to) 

*  —  oo 

X  £(00)[(xi  +  *2)/2.  0>]e~iuTdw.  (D4) 

The  signals  received  by  the  antennas  are  filtered  so  that  only  those  radio  waves  with  frequency  to 
within  a  narrow  bandwidth  Ato  centered  on  a  median  frequency  contribute  to  the  data.  Thus  upon 
transforming  Eq.  (D4)  by  setting  to  =  to'  +  to  we  have 

—  Aw  A 

r(0>(x|,  x2,  r)  =  e -'"T/(2x)3  J  G(oo)(x{  -  x*.  to'  +  5) 

x  £<0D)[(x,'  +  x2)/2,  to'  +  u]  [(to'  +  5)/c]e"'"'Vto'.  (D5) 

Next  we  make  the  quasi-monochromatic  approximation 

Ato  «  to,  (D6) 
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so  that  u'  «  u  over  the  range  of  integration  in  Eq.  (D5).  We  can  set  (w'  +  «)/c  =  w/c  =  k 
where  the  left-hand  side  of  this  expression  does  not  appear  in  an  exponent  in  Eq.  (D5).  Now  con¬ 
sider 


r 


G(0D)(xj  +  Xj,  u'  +  w)  = 


||  (Rs)/cos2  6}  exp  li[(u'  +  u>)/c]  s  •  (xj  -  x£))d2 s.  (D7) 


If  we  restrict  the  size  of  the  antenna  so  that 

|  xj  Xj  |  «  (2x c/A  u),  (D8) 

then  we  may  neglect  the  effect  of  u'  in  Eq.  (D7)  so  that,  if  we  assume  that  /(00)(x)  is  independent  of 
w',  and  we  can  write  Eq.  (D5)  as 

r<0)(xj,  xj,  r)  =  (*/2t)2G<od)(x;  -  Xj,  5»<->[(*i  +  x})/2,  r]  exp  (-i«r),  (D9) 


where 


^<00)(x;,r)  4  (1/2t)  I*  ^(xi.S + 


(DIO) 


To  eliminate  the  effects  of  partial  coherence  of  the  source  on  the  data  we  normalize  the  mutual  coher¬ 
ence  function  in  the  usual  way  to  obtain  the  complex  degree  of  spatial  coherence 


7<0)(xj,  «5,  T)  Ar(0)(xj,  X*.  T)/Vr(0)(xj,  xj,  r)r(0)(xi,  x;,  T) 

=  G(oo)(xj  -  x*.  <3)m  "’[(xj  +  x$)/2,  T]/  [c(oe)(0,  5)  V^'Vxj,  r) 


(DU) 


Since  /(0O)  (as  well  as  /?[,  R2  and  6)  is  a  slow  function  relative  to  then  it  follows  from  the  prop¬ 
erties  of  Fourier  transform  pairs  that  Ji  is  a  slow  function  relative  to  Aro°\  The  temporal  Fourier 
transform  in  Eq.  (DIO)  does  not  change  this  relationship  so  that  we  can  assume  that 
G(00)(xj  -  x{)  =  0  except  over  very  small  domains  in  which 

r’uxj  +  *0/2,  T)  =  rw  T)  ~  A(W)(X2,  r).  (D12) 

Upon  substituting  from  Eq.  (D12)  into  Eq.  (Dll)  we  have 

7<°>(xj,  x£,  r)  =  7<0)(xj,  Xj,  0)  =  G(oe)(xj  -  x*.  «)/G(”>(0,  5),  (D13) 

so  that  upon  substituting  from  Eq.  (D2),  and  inverting  the  Fourier  transform  Eq.  (D13)  becomes 


l(Rs)  =  G0(cos2  6/\2R2)  |J  g(xL,  o>)  exp  (-fits  •  x'_)d2x'_,  (D14) 
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where  G0  =  G(ao)( 0,  w)  is  a  constant  and 
g  (*i  ~  *2.  «)  ^  7®(*{.  *2.  0), 
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iffi  |  xj  |  <  a  and  |  x£  |  :£  a , 


=  0  otherwise,  (D15) 

is  a  stationary  correlation  coefficient.  Equation  (55)  follows  by  comparing  Eq.  (D14)  to  Eq.  (52)  by 
using  Eqs.  (Dll)  and  (D15).  Since  the  function  can  be  measured  only  over  the  portion  of  the 
z  —  0  plane  that  is  inside  of  the  circle  of  radius  a  it  is  necessary  to  define  g  to  be  zero  for 
I  *i  ~  *2 1  >  a  t°  invert  the  Fourier  transform  in  Eq.  (D13)  to  obtain  Eq.  (D14). 
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CALCULATION  LEADING  TO  EQUATION  (60) 

The  derivation  of  Eq.  (60)  is  virtually  identical  to  that  of  Eq.  (55)  except  for  one  slightly  dif¬ 
ferent  assumption  and  some  different  coordinates  and  factors.  The  cross-spectral  density  function 
over  the  antenna  hemisphere  in  Fig.  1 1  is  given  by 

H'i-’fl?,*.  /? 2*2)  =  (*/2t)2(cos  6\  cos  B2/R\R2 )  exp  [ik  (R  ,  -  R2) 


where 


X  7(S|  —  S2,  w)  ft  [(S[  +  *2)/2,  «], 


7<0)(s_,  w)  =  JJ  /(0)(x',  o>)  exp  ( — /A s _  •  x'^hi' , 


(El) 


(E2) 


and 


£(0)(s+, «) 


00 

J|  /i<0)(x_,  «)  exp  (~iks+  ■  x_)d1x_ 

-  Oo 


(E3) 


(see  Section  2,  Eqs.  (33),  (34),  and  (35)),  where  /<0)  is  the  intensity  distribution  over  the  z  =  0 
source  plane,  and  n(0>  is  the  complex  degree  of  spectral  coherence  over  the  source  plane.  Equations 
(El)  through  (E3)  are  different  from  Eqs.  (Dl)  through  (D3),  but  this  does  not  affect  the  derivation  of 
Eq.  (60)  in  general. 

Again  as  in  Appendix  D  we  assume  that  the  source  is  quasi-homogeneous  so  that  /(0)  is  a  slow 
function  relative  to  nl0).  The  derivation  then  proceeds  just  as  in  Appendix  D,  with  the  same  assump¬ 
tions,  until  we  come  to  the  requirement  that  the  antenna  plane  be  restricted  so  that  Eq.  (D8)  is  satis¬ 
fied.  Instead  of  the  assumption  given  by  Eq.  (D8)  we  must  make  the  similar  assumption  that  7(0)  van¬ 
ishes  outside  of  the  domain  where 


I*' I  <  <  Or c/Aw).  (E8) 

With  this  assumption,  and  with  the  assumption  that  the  source  intensity  /(0)  is  effectively  constant  over 
the  frequency  domain  passed  by  the  receivers,  we  may  remove  /<0)(s,  w)  from  the  integration  over  01' 
as  we  removed  G(a>)  from  a  similar  integral  in  Appendix  D. 

The  derivation  then  proceeds  as  before  until  we  finally  obtain  instead  of  Eq.  (D14) 

OO 

7((V,  w)  =  70/X2  ||  g'(s_)  exp  (i*s_  x')d2 s_,  (E14) 
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where,  instead  of  Eq.  (D15), 

g'( St  -  §2)  4  r<->(tflSl,  /f2S2,  0)  exp  -  ^2)]/Vr(*)(/f1s1,  Rysu  0)r(*)(/?2s2,  /?2S2,  0), 

iffi  SJ-,  +  s2y  <  1  and  sh  +  s2y  s  1, 

=  0  otherwise,  (H15) 


is  a  stationary  correlation  coefficient,  and  where  IQ  is  a  constant  given  by 


/0  =  ff  /(£V,  w)dV. 


(E16) 


By  comparison  of  Eqs.  (E14)  and  (E15)  with  Eqs.  (58)  and  (59)  it  is  clear  that  they  are  equivalent  if 
g'  is  given  by  Eq.  (60).  Because  g'  is  defined  physically  only  for  unit  vectors  s,  and  S2  pointing 
from  the  origin  toward  the  z  >  0  half  space  so  that  j1x  +  sfy  <  1,  and  s£x  +  s2y  rS  1,  it  is  neces¬ 
sary  to  define  g'  =  0  for  all  other  values  of  su,  and  s2jt ,  s2y  as  indicated  in  Eq.  (E15)  to  invert 
a  Fourier  transform  to  obtain  Eq.  (E14). 
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